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Abstract 

Wc present the first measurement of clustering properties of low mass galaxies with a stellar mass down 
to M* ~ 10 9 M o at 1 < z < 4 in 24.4 arcmin 2 of the GOODS-North region with a depth of K AB ~ 25, based 
on the near infrared observations performed with MOIRCS at the Subaru Telescope. The correlation 
amplitude strongly depends on the K-band flux, color, and stellar mass of the galaxies. We find that 
A'-band luminous galaxies have a larger correlation length than A-band faint galaxies. For color selected 
samples at 2 < z < 4, distant red galaxies with J — K > 1.3 show a large bias of b ~ 7.2 ± 1.3 on scales up 
to 9 ~ 100" or 3.1 comoving Mpc, while blue galaxies with 0.5 < J — A < 1.3, in which most Lyman break 
galaxies are populated, have a weak clustering signal on large scales, but a possible strong small scale 
excess at 8 < 10". For massive galaxies with M* > 10 10 M Q , we estimate a correlation length and bias to 
be ro ~ 4.5 h _1 Mpc and b = 1.9 — 3.5, which are much larger than those of low mass (M* ~ 10 9 -10 10 Mq) 
galaxies. The comparison of our measurements with analytic CDM models constrains the properties of 
hosting dark halos, and indicates that the low mass galaxies would be progenitors of galaxies with a typical 
luminosity of L < L* in the local Universe. The blue galaxies in low mass samples are more strongly 
clustered in more massive halos with higher occupation numbers than low mass red galaxies. This fact 
suggests an environment effect due to the halo mass on star formation activity at high-z. 

Key words: cosmology: observations — galaxies: evolution — galaxies: high-redshift — infrared: 
galaxies — cosmology: large-scale structure of universe 



1. Introduction 

The spatial correlation of galaxy distributions with that 
of the underlying dark matter is one of fundamental in- 
formation to understand galaxy evolution. Under a fixed 
cosmological framework of cold dark matter models, the 
comparison of observed clustering properties of galaxies to 
the theoretical predictions allows us to obtain, in a statis- 
tical manner, typical masses of the dark matter halos host- 
ing galaxies even at the high-z Universe (e.g., Moustakas 
& Somerville 2002; Ouchi et al. 2004; Hamana et al. 2006). 
Therefore the evolution of galaxy clustering strength and 
stellar mass assembly in galaxies as a function of redshift 
will give us clues to the assembly history of not only the 
stellar mass but also dark halos. Thus, it allows us to 
study ancestor-descendant connections of galaxies at dif- 
ferent redshifts. 

Previous studies of galaxy evolution at high-z have 



been limited mainly to massive galaxies with stellar mass 
M* > 10 10 M Q . According to hierarchical clustering mod- 
els, low mass galaxies played an important role in the 
high-z Universe as building blocks, merging into massive 
galaxies. Therefore the comparison of the clustering prop- 
erties of massive and lower mass galaxies will give us a 
more general view of galaxy evolution in dark matter. In 
this context, we have performed deep near-infrared (NIR) 
imaging observations (Kajisawa et al. 2006) to measure 
the properties of high-z galaxies based on a A-selected 
catalog, focusing on low luminous (or low mass) galaxies 
in comparison with massive galaxies. 

Clustering properties are well studied at z < 1 in optical 
band with spectroscopic redshift data (e.g., Zehavi et al. 
2002, 2005; Norbcrg et al. 2002; Coil et al. 2004; Meneux 
et al. 2006; Li et al. 2006). For higher redshift, where 
spectroscopic redshift data arc limited to bright galaxies, 
multicolor selection techniques have efficiently revealed an 
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abundant population of galaxies. Optically selected UV 
luminous Lyman break galaxies (LBGs) (Steidel et al. 
1996) and NIR selected galaxies, such like distant red 
galaxies (DRGs) (Franx et al. 2003) are among the galax- 
ies selected by such color selection techniques for high-z 
galaxies. It is important to notice that the clustering prop- 
erties of galaxy samples would be sensitive to the selection 
criteria. LBGs are generally biased to UV luminous ac- 
tive star-forming galaxies, and may give biased samples in 
terms of stellar mass (van Dokkum et al. 2006). On the 
other hand, galaxy catalogs compiled with NIR data allow 
us to construct a mass-selected sample because NIR emis- 
sion is less affected by dust extinction and closely traces 
total stellar mass. Thus a NIR selection may give a suit- 
able sample to study the clustering properties of galax- 
ies based on the stellar mass. Nonetheless, the previous 
studies at high-z on the basis of NIR observations have 
been limited in deep (K < 26), but small field of view (4.5 
arcmin 2 ) (Labbe et al. 2003) or in shallow (K < 23.5) and 
wide fields (> 30 arcmin 2 ) (e.g., Glazebrook et al. 2004; 
Quadri et al. 2007; Foucaud et al. 2007; Grazian et al. 
2006b). 

To examine the clustering properties of galaxies at 
1 < z < 4, we will analyze in this study the medium 
deep (K < 25) and wide-field (28 arcmin ) NIR data 
taken with MOIRCS on the Subaru Telescope (Kajisawa 
et al. 2006) and a publicly available data in the Great 
Observatory Origins Deep Survey North (GOODS-N) re- 
gion (Giavalisco et al. 2004). Although the depth (K = 
25.1 at the 90% completeness) is shallower than K = 25.7 
of Daddi et al. (2003), the field of view is 5.5 times larger. 
The limiting magnitude, ~1.5 mag deeper than that of 
e.g., Grazian et al. (2006b), will give us data statistically 
robust in number and with smaller field variance for the 
clustering properties for low luminous (or low mass) galax- 
ies. Applying analytic models for the spatial clustering of 
dark matter to the observational results, we estimate the 
mass of dark matter halos hosting massive and low mass 
galaxies from the bias of galaxy-dark matter distribution 
(Bullock et al. 2002; Moustakas, Somerville 2002; Allen 
et al. 2005). 

The present paper is organized as follows. In §2, we 
present a brief account of the observations and data re- 
duction. The photometric redshift and stellar mass of 
observed galaxies are obtained. We describe in §3 the an- 
gular clustering estimates, and the dependence on color, 
flux, and mass of galaxies at z = 1-4. Summarizing the 
result in comparison with previous studies, we estimate in 
§4 the galaxy-dark matter biases to investigate the proper- 
ties of dark halos hosting high-z galaxies and the relation 
between galaxies and dark halos. In §5, we discuss our 
results and relevant literature in conjunction with galaxy 
evolution scenarios and possible candidates of the descen- 
dant in the local Universe. 

Throughout this paper, we assume O m = 0.3, Q\ = 0.7, 
and Ho = 70 km s -1 Mpc -1 . To facilitate the compar- 
ison to previous studies, the results of the correlation 
length (ro) are expressed using h = 1 (Hq = 100 km s _1 
Mpc -1 ). We use the AB magnitude system (Oke, Gunn 



1983; Fukugita et al. 1996), unless otherwise stated. For 
comparison with previous observations in literature, we 
convert K in Vega system using Kab = -KVega + 1-82. 

2. Data 

2.1. Observations and Data Analysis 

We performed MOIRCS Deep Survey (MODS), J, H, 
if s -band imaging observations in GOODS-N region with 
the Multi-Object InfraRcd Camera and Spectrograph 
(MOIRCS; Ichikawa et al. 2006) on the Subaru Telescope 
in April and May, 2006. MOIRCS has a field of view 
of 4x 7 arcmin 2 with 0".117 pixel scale. We observed 
one field of view of MOIRCS centered at 12 h 36 m 46?62, 
+62°13T5."6 (J2000), which includes the Hubble Deep 
Field North (HDF-N, Williams et al. 1996) at the cen- 
ter. The sky was clear and stable during the observations 
under good seeing condition. We choose in this study 
highest-quality data sets which include only the frames 
with seeing size smaller than 0".5 (FWHM). 

The data were reduced in a standard manner using the 
IRAF software package. The details of the data processing 
and the data quality will be described at length in Konishi 
et al. (in preparation). FS23 and FS27 in the UKIRT 
faint standard stars were used for the flux calibration. 
The FWHMs of the point spread function (PSF) for the 
combined images are 0".42, 0".41, and 0".40 in J, H, K s - 
bands, respectively. The total exposure times of the final 
images are 5.0 (J), 1.8 (H), and 7.7 (K s ) hours. 

The source detection is performed on the i^ s -band im- 
age using the SExtractor image analysis package (Bertin, 
Arnouts 1996). We use MAG.AUTO of the SExtractor as 
the total K magnitudes of the detected objects. In total, 
we have detected 2564 objects down to if ~ 26.3. The 
5a limiting magnitudes in 0".85 aperture arc J = 25.7, 
H = 24.2, and K = 25.3. 

The photometric errors and the completeness of 
the detection in K s band are obtained by Monte 
Carlo simulations with mock images produced by the 
IRAF/ ARTD ATA package. We mingle many simulated 
galaxies in our image and then search them using the same 
photometric method used to detect actual galaxies. We 
find the 90% completeness at K = 25.1. We also test sen- 
sitivity to false detections by running SExtractor on the 
inverted _fT s -band image. Only 15 spurious objects are ex- 
tracted at K < 25, which is reasonably negligible for the 
present study. 

The field of view of MOIRCS is divided into two 4' x 3'.5 
fields and focuses on two separate focal plane arrays 
(Ichikawa et al. 2006). Excluding the edge area of each 
frame, where the S/N is lower than the center due to 
dithering or pixel defect, we obtain the images with the 
K ~ 25.0 limiting magnitude uniform over the field of 24.4 
arcmin 2 , which is smaller than the original observation 
area (28 arcmin 2 ). The final catalog used for the present 
study contains 1959 galaxies to the 90%completeness limit 
K = 25. 

Our survey area is covered by deep ACS images with 
the F435W, F606W, F775W, and F850LP bands of HST 
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circles, open circles, and crosses, respectively. The typical 
sp error for each z bin is shown by a thin line. 



Fig. 1. Upper panel: the spectroscopic (z sp ) vs. photometric 
( 2 photo) redshifts for 527 galaxies in the present region. Lower 
panel: relative scatter (z photo — z sp ) / ( 1 + z sp ) as a function of 
z Bp . The inset histogram shows the distribution of the photo- 
metric redshift error. The error distribution is approximated 
by a Gaussian centered at 0.004 with an rms of 0.06 (solid 
line). 
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Fig. 2. Redshift distributions for 1959 galaxies with 
K < 25.0. The spectroscopic samples for 527 galaxies are 
depicted by the filled histogram. 



in the GOODS-N region (Giavalisco et al. 2004). (For 
convenience, we designate them as B, V, i, and z when 
referring to the filters, respectively.) A deep £/-band im- 
age taken with KPNO 4m telescope is available in Capak 
et al. (2004), though the PSF is much larger than those 
of ACS and MOIRCS. Using our MOIRCS data and the 
public dataset, we have produced a high-quality multi- 
color catalog of galaxies in the MODS region. 

The ACS images and our K s and H images are con- 
volved with a Gaussian kernel to match the PSF to 0".42 
of J band, which has the largest PSF in the present image 
set, except U band. For the color measurements, we use 
0".85 aperture (2xFWHM). For the [/-band image with 
FWMH=1".3, we measure the flux in 2".8 aperture. The 
aperture correction is made using the difference of the pho- 
tometric measurements with apertures 0".85 and 2". 8 on 
the B-band image after convolved with a Gaussian kernel 
to match the PSF to 0".42 (J) and 1".3 (U), respectively. 

2.2. Photometric Redshift and Stellar Mass 

We compile spectroscopic redshifts for 527 galaxies in 
the present region from the literature (Wirth et al. 2004; 
Cowie et al. 2004; Barger et al. 2003; Cohen et al. 2000; 
Cohen 2001; Treu et al. 2005; Erb et al. 2004; Steidel et al. 
2003; Dawson et al. 2001; Reddy et al. 2006), excluding 
several unreliable identifications. (Note that the spectro- 
scopic redshifts arc available for only 7 DRGs among 115 
DRGs in the present region.) For all the detected ob- 
jects, we estimate their redshifts photometrically with the 
photometric data of UBVizJHK by fitting the observed 
flux with model SEDs, following the method described in 
Kajisawa, Yamada (2005). We use the model SEDs of 
Bruzual, Chariot synthetic library (GALAXEV; Bruzual, 
Chariot 2003), the Calzetti extinction law (Calzetti et al. 
2000), and Hi absorption (Madau 1995). Free parameters 
used in the fitting are redshift, spectral type, age, and ex- 
tinction. We assume a star formation rate (SFR) that 
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Fig. 4. Sky distributions of the galaxies at (a) 1 < z < 2 and 
(b) 2<2<4. The galaxies with J-K> 1.3, 0.5 < J -if < 1.3 
and J — K < 0.5 are depicted with filled circles, open cir- 
cles, and crosses, respectively. The large and small symbols 
represent the galaxies with stellar mass M* > 10 10 M o and 
M* = 10 9_10 M Q , respectively. 

decays exponentially with time, SFR oc e~*/ T , where r is 
the time scale between 0.01 and 30 Gyr. The metallicity 
is changed from 0.005 to 2.5 solar values. The initial mass 
function (IMF) of Chabrier (2003) is adopted with upper 
and lower mass cutoffs of mi = 0.1M Q and m u = 100 AIq. 

We check the accuracy of our photometric redshift mea- 
surement using the objects with available spectroscopic 
redshifts. Figure 1 compares the spectroscopic and photo- 
metric rcdshifts. There are several outliers at z sp ~ 0.5 and 
2 < z < 4, for which difference is as large as \Az\ ~ 2. The 
error distribution (the inset of the lower panel of Fig. 1) is 
nearly Gaussian and scatters around zero with an rms er- 
ror of Az/(1 + z sp ) ~ 0.06. Excluding \Az\/(l + z sp ) > 0.5, 
we obtain the rms error ~ 0.12, which is comparable with 
the results 0.12, 0.09, and 0.05 of Quadri et al. (2007), 
Rudnick et al. (2003) and Grazian et al. (2006a). The 
outliers are possibly caused by the confusion of Lyman 
break and 4000A/Balmer break in the photometric red- 
shift technique. We will discuss the influence of the out- 
liers to the evaluation of clustering properties of galax- 



ies in §3. Figure 2 shows the redshift distribution of the 
present if-band selected galaxies. We hereafter use the 
spectroscopic redshifts whenever available. 

GALAXEV gives the stellar mass-to-light ratio (M»/L) 
and rest frame colors for each template, so that we can ob- 
tain the total stellar mass (M») of the galaxies using the 
total absolute magnitude in rest frame V band, corrected 
for dust extinction. Note that our stellar mass estimate 
with the Chabrier (2003) IMF is systematically about 1.8 
times smaller than that with the Salpetcr IMF (Salpeter 
1955). We plot the mass distributions as a function of 
redshift in Fig. 3, where the galaxies are divided into 
three groups with colors J — K > 1.3, 0.5 < J — K < 1.3, 
and J-K < 0.5. All of 52 LBGs of Steidel et al. (2003), 
which are located in the present region, are identified in 
our K-selected catalog. The colors of most LBGs are dis- 
tributed in 0.5 < J — K < 1.3. Therefore, we divide the 
bluer galaxies into two groups by J — K = 0.5 for the fol- 
lowing analysis. We see that most DRGs (J — K > 1.3) 
are massive galaxies (3.6 x 10 10 Mq on average) residing 
at 2 < z < 4. 

The sky distributions of the sample galaxies are repre- 
sented in Fig. 4, where we see, as a visual impression, 
the difference of the distributions between different sam- 
pling criteria by the redshift, color, and stellar mass. In 
what follows, we examine the differences in a quantitative 
manner. 

3. Angular Correlation Function 

We quantitatively measure the clustering properties of 
galaxy distributions using the angular two-point correla- 
tion function (ACF) w{9) (Peebles 1980). To examine 
the dependence on galaxy properties, we define a variety 
of different subsamplcs based on the redshift, flux, color, 
and stellar mass (see the following section for the defi- 
nitions). We adopt the minimum variance Landy-Szalay 
estimator (Landy, Szalay 1993): 

RR(9) 



t ^_ DD{9)-2DR{9) 
W{d) = RkW 



(1) 



where DD(9) is the observed number of galaxy pairs with 
separation 8, DR{9) the number of pairs between the ob- 
served galaxies and random samples, and RR{9) the num- 
ber of pairs in the random catalog. In the calculation, 
we distribute the same number of random samples as of 
the observed sample in the same geometrical constraint 
of the observations. To reduce the noise in the random 
pair counts, we repeat the measurement until the random 
samples count over 100,000. We compute w(9) in log- 
arithmic bins of width Alog(9) = 0.4. The ACF errors 
are estimated by a bootstrap re-sampling technique (e.g., 
Ling et al. 1986). Note that our error estimate does not 
include field variance. 

The ACFs measured in a finite sky region are influ- 
enced by the, so called, integral constraint IC, which is 
expressed by 



IC' 



1 



dn 1 dn 2 w T (9), 
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Fig. 5. Upper panel: The ACFs, w(8), of the if-band selected galaxies at (a) 1 < z < 2 and (b) 2 < z < 4. The filled circles show the 
observations with la bootstrap error. The solid curve with dotted la error range is the best-fit power law AO-P with fixed /3 = 0.8 
for 1"— 100", corrected for integral constraint IC. The dash curve shows the ACF of dark matter predicted by the nonlinear model 
of Peacock, Dodds (1996). Lower panel: Galaxy-dark matter bias b, defined as b{6) = *J w(6)/wt}m ■ The thin line with dotted lines 
of ±lcr error indicates the average bias and the dispersion in the bins with 8" < 8 < 100". 



where wt(9) is the true ACF and f2 is the solid angle of the 
field (Groth, Peebles 1997). We assume a power-law ACF 
of the form, w(9) = A9~ 13 , and determine its amplitude, 
A, by fitting the function, 

w{6) = AO' - IC, (3) 

to the observations. Our measurement, however, cannot 
significantly constrain the slope of the correlation function 
because of the relatively small number of pairs and small 
angular range. Therefore, we take a fixed slope of f3= 0.8. 
The amplitude A is computed by fitting the ACF over the 
range 1 — 100" because the finite size effect is serious for 
larger separation and our aperture photometry size 0".85 
would hardly separate close pairs in 9 < 1", if any. We 
estimate the integral constraint numerically as 



IC 



(4) 



following Roche ct al. (2002). 

Since ACF quantifies the clustering properties of galax- 
ies projected on the sky, it reflects the combination of the 
rcdshift distribution of selected galaxies and the cluster- 
ing in three-dimensional space. The spatial correlation is 
usually expressed by a power law as 



€ = (r/ro)" 



(5) 



where r is the spatial separation between objects, tq the 
correlation length, and 7 the slope of power law. The 
slope (3 of ACF is related to 7 with 7 = /? + 1. We use the 
Limber equation (Limber 1953) to infer tq from w(9) in 
comoving units: 



w{9) 



(6) 
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Fig. 6. The dependence of the clustering length on K mag- 
nitude at 1 < z < 2 (left) and 2 < z < 4 (right). Open circles 
represent the results of the present study. Square and triangle 
show the results of Daddi et al. (2003) (FIRES) and Quadri 
et al. (2007) (Q07), respectively. 
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where A Q bj is the number of galaxies in a sample and 
dN / dz is the rcdshift selection function of sampled galax- 
ies. For dN/dz, we use the actual redshift distribution of 
each sub-sample. 

Since the scatter in the photometric redshifts is very 
small as shown in the previous section, the error in the 
photometric redshifts does not have a significant influence 
on our results. The outliers of the photometric rcdshift act 
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Fig. 8. The results of the correlation length, ro, for red 
galaxies based on X-selected samples at z > 1. The filled cir- 
cle represents our result for DRGs ( J — K > 1.3) at 2 < z < 4. 
The other symbols show the previous studies, as labeled in 
the legend: Quadri et al. (2007) (Q07), Foucaud et al. (2007) 
(F07), and Grazian et al. (2006b) (G06) for DRGs, Daddi 
et al. (2003) for J - K > 0.7 galaxies (FIRES), and Daddi 
et al. (2002) for ERO. 

as contaminants. If the contaminants have a uniform dis- 
tribution, the measured amplitude A should be multiplied 
by 1/(1 — f c ) 2 , where f c is the fraction of the contami- 
nants. As shown in the previous section, the outliers at 
z < 1 could be counted at the redshift range 1 < z < 4 
due to the photo- z error. The numbers of the interlopers 
observed at 1 < z < 2 and 2 < z < 4, which are actually 
located at z < 1, arc ~1% and ^3% of all galaxies, re- 
spectively (Fig. 1). Among the 877 galaxies at z < 1, 358 
have reliable spectroscopic redshift. Therefore, out of the 
other 498 galaxies, ~5 and ~15 galaxies are potential in- 
terlopers which could be counted in the redshift range of 
1 < z < 2 and 2 < z < 4, respectively. They are 1% and 3% 
of the samples. Therefore, the photo- z error decreases the 
amplitude by ~6% at most. It corresponds to an error of 
~ 4% on ro , which is much smaller than the fitting error 
of r . 

We also estimate the bias parameter &, which is defined 
as the ratio of the clustering amplitudes of galaxies to that 
of dark matter, 



w= Mo) 



WDM 



(10) 



where w{&) and u>dm are the observed ACF corrected for 
IC and that of dark matter computed by the Limber's 
projection with the nonlinear model of Peacock, Dodds 
(1996). We calculate wdm using the same observational 
selection function of each subsample. Note that w(9) of 
Eq. 1 can be negative due to a statistical error. In such 
cases we artificially set b(9) = ~\/\w(8)\/wbm ■ 

4. Results from Clustering Measurements 

We present results of clustering measurements for sub- 
samples in the present ^-selected galaxy catalog with 
K < 25. We define subsamples by (1) K flux, (2) J — K 



color, (3) stellar mass M* and (4) rest frame U — V color. 
Also, in the ACF analysis, we divide the catalog into two 
redshift bins of 1 < z < 2 and 2 < z < 4 to investigate the 
redshift evolution of the galaxy clustering. The redshift 
intervals are chosen so that each sample has roughly the 
same number of galaxies (~ 500) and similar time intervals 
(~2Gyr). 

In Fig. 5, the observed w(B) and the best fit power-law 
model with a fixed slope of j3 = 0.8 and galaxy-dark mat- 
ter bias b are depicted. It is found from the figure that 
the clustering bias of the present sample increases with z. 
In the higher redshift sample, the bias excess at the small 
scales (9 < 8") is distinctly observed. In the galaxy catalog 
with flux limit, various populations of galaxies are likely 
to be mingled. Populations uniformly distributed in space 
would smear out actual strong clustering of some popu- 
lations, if any. To see what populations are contributed 
to the large bias or small scale excess, we disentangle the 
populations using different selection criteria. 

4-1- Clustering on K-band Flux Selection 

We show the dependence of the clustering length on K 
magnitude in redshift bins of 1 < z < 2 and 2 < z < 4 in 
Fig. 6, where the ACF analysis is made with the limiting 
magnitudes of K = 23, 24, 25 for our sample. In this figure, 
we see that the clustering length increases with K-band 
luminosity in both redshift bins; more luminous galaxies 
tend to be more strongly clustered. The result conflicts 
with the previous study by Daddi et al. (2003) based on 
the K ~ 26 sample in 4.5 arcmin 2 HDF-S field, in which 
they concluded that the clustering length declined only 
slightly from K ~ 21 and remained as high as ro ~ 5 down 
to K ~ 26. Fluctuations in the clustering amplitude due 
to field variance could affect their estimate derived in the 
small HDF-S field. 

Quadri et al. (2007) studied in 300 arcmin 2 the angular 
correlation functions of K selected galaxies with 2 < z < 
3.5 and K < 22.8. Their clustering measurement based 
on the shallow but wide-field data is statistically more 
robust than ours at the bright magnitude of K < 23. The 
correlation length for their sample is ro = 6.0^°'^ h~ 1 Mpc 
with 7 = 1.8, which is in good agreement with our result 
ro = 4.8^48 ^ _1 Mpc with the same selection criteria (K < 
23), even though our error is much larger due to the small 
area and small number (N = 54) of galaxies. The bias 
b = 4.1 ± 2.0 for our K < 23 sample is also consistent with 
that of Quadri et al. (2007) (b = 3.3 ±0.5) within la error. 
They found no dependence of the correlation lengths on K 
magnitudes, whereas our deeper observation have revealed 
the strong dependence on K flux at a fainter magnitude of 
K > 23. The result that at high-z the dependence of the 
correlation lengths on K flux is weak for bright galaxies 
(Quadri et al. 2007), whereas that for faint galaxies is 
strong (this study), suggests opposite tendency for the 
galaxies at the local Universe (e.g., Norberg et al. 2002; 
Li et al. 2006) or for LBGs at z - 4 (Ouchi et al. 2005). 
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Fig. 7. Same as Fig. 5, but for DRGs (J — K > 1.3) and bluer galaxies at 2 < z < 4. 



Clustering of DRGs and Bluer Galaxies at 2 < z < 4 

To study the color dependence of the clustering prop- 
erties, we present the results of ACF analysis based on 
J — K color selection in Fig. 7. It is found from Figure 7 
that the bias increase with the J — K color; the stronger 
clustering for the redder J—K color. 

The clustering lengths of DRGs are measured by the 
previous studies (Quadri et al. 2007, Foucaud et al. 
2007, Grazian et al. 2006b). We compare their re- 
sults with ours in Fig. 8. As shown in the figure, 



our result (ro = 9.2 
J - K > 1.3 at 
within lcr error (ro = 
2007; r = 11.1±|-| hr 1 

i+3.0 



Jj'jjj h J Mpc for the sample with 
2 < z < 4) agrees with their results 
12.0±?;g /i _1 Mpc, Quadri et al. 
Mpc, Foucaud et al. 2007; r = 



13.4+g;2 /i _1 Mpc, Grazian et al. 2006b). We note, how- 
ever, that a marginal tendency of a larger correlation 
length for a brighter sample is observed among the re- 
sults. 

Daddi et al. (2003) found the color segregation of the 
clustering at 2 < z < 4, dividing the catalog into J—K > 1.7 
and bluer galaxies so that the number of the galaxies 
become similar. (J — K = 1.7 in Vega corresponds to 
J - K = 0.74 in AB of MOIRCS photometric system.) 
Although their sample is 0.7 mag deeper than our data 
at the same 90% completeness level, our catalog is ex- 
pected to well sample the galaxies with J — K > 0.74 at 
2 < z < 4. In fact, the surface number density 11.0 ±0.7 
galaxies arcmin -2 in our catalog is very consistent with 
their sample (10.9 ±1.6 arcmin -2 ). To check the con- 
sistency with the result of Daddi et al. (2003), we sub- 
sample the galaxies of our catalog by J — K = 0.74 and 
measure the correlation lengths. For the redder sample 
(J-K > 0.74), we obtain r = 5.3±° l /i _1 Mpc, which 
is much small than the value provided by Daddi et al. 
(2003) (8.3 ± 1.2 /i -1 Mpc). For the bluer sample, we ob- 
tain ro = 2.h + } 2 \ /i _1 Mpc, which is slightly smaller than 
that of Daddi et al. (3.5±3 £ h _1 Mpc), though the dif- 
ference is within 1 a error. The discrepancy may likely 
arise from the ficld-to-field variation as discussed in the 



previous subsection. 

In addition to the J — K color segregation of the clus- 
tering, it is important to notice that the sample with 0.5 < 
J— K < 1.3, which include the most of LBGs, clearly shows 
the excess clustering at small scales (9 < 8"), while the 
sample has a weaker clustering amplitude at 8" < 9 < 100". 
Interestingly, the correlation length (ro = 4.lti'g/i _1 Mpc) 
and bias (b = 2.6 ± 0.4) for 0.5 < J - K < 1.3 galaxies 
are consistent with those of LBGs for L> L* at (z) = 3 
{r =5.0±l 7 7 h- 1 Mpc, 6 = 2.7±0.4) (Giavalisco, Dickinson 
2001; Ouchi et al. 2004). 

4-3. Clustering on Stellar Mass Selection 

Next, we examine the stellar mass dependence on the 
galaxy clustering. We divide the sample into the two mass 
ranges of massive (M* > 10 10 M@) and low mass (M* = 
10 9-10 Mq) populations. The result is presented in Fig. 
9, which shows that massive galaxies are larger biases than 
the low mass galaxies at the same redshift ranges. 

We examine which populations mainly contribute to the 
massive galaxies. Some massive galaxies at 2 < z < 4 would 
be overlooked by the standard DRG selection technique 
as shown in Fig. 3 (see also Conselice et al. 2007 and 
van Dokkum et al. 2006). We count 124 galaxies with 
M* > 10 10 M Q at 2 < z < 4. Among them, 54 galaxies 
are found to be DRGs with J — K > 1.3. van Dokkum 
showed that the fraction of DRGs at 2 < z < 3 is 69% by 
number among the massive galaxies with M* > 10 11 M Q , 
using the Salpeter IMF. At the same 2 < z < 3, we count 
16 galaxies with M* > 6 x 10 10 Mq which corresponds 
M* > 10 n M Q with the Salpeter IMF. Among them, 13 
galaxies (81%) arc DRGs; the fraction is consistent with 
the result of van Dokkum et al. (2006), if the field variance 
is taken into account. LBGs are another population of 
massive galaxies (e.g., Rigopoulou et al. 2006). In fact, 
van Dokkum et al. (2006) studied the population of LBGs 
with A/* > lO n M at 2 < z < 3. However, it is noted 
that the 52 LBGs of Steidel et al. (2003) in the present 
region have various stellar masses with the wide range of 
7 x 10 8 M Q - 9 x 10 10 M o with no typical mass. Among 
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them, 12 LBGs have masses of M* > 10 10 M Q . 

Counting DRGs and LBGs, we find that only 53% of 
massive galaxies are classified by the combination of the 
color selection techniques. In summary, massive galaxies, 
including the DRGs sample (J — K > 1.3) and a minor 
portion of LBGs in 2 < z < 4, contribute large bias at 
8 > 8" (or > 0.25 Mpc at z ~3) with large correlation 
length. The redder massive galaxies are more strongly 
clustered as discussed below. 

4-4- Clustering on the rest frame U—V color 

Finally, we examine the rest frame U — V color depen- 
dence on the galaxy clustering. We divide the sample into 
the two subsamples by the rest frame U — V = 0.7. Note 
the color U — V ~ 0.7 in AB magnitude corresponds to 
the color of AO type stars (U — V ^ in Vega). Therefore, 
U — V < 0.7 galaxies are likely to be in active star forma- 
tion. Also, it should be noted that most of the massive 
galaxies (M* > 10 10 M©) at 1 < z < 2 have color redder 
than U —V = 0.7. 

We summarize the results in Table 1. The first column 
is the redshift range, in which we examine the cluster- 
ing properties for each subsample. The second column 
is selection criteria for subsampling. The number of the 
galaxies, N, for each subsample is shown in the third col- 
umn. The average redshift (z) and stellar mass (M*) arc 
calculated for each subsample in column 4 and 5, respec- 
tively. The sixth column is the best fit parameter A in Eq. 
3 and the la error with a constrained slope (3 = 0.8. Using 
the Limber equation (Eq. 7), we convert the parameter to 
comoving correlation length, ro, in units of ft -1 Mpc in the 
seventh column. The eighth column is the bias averaged 
in the bins, where we use the bins with 8" < 8 < 100" to 
avoid a small scale excess. In last column, we compute 
the number density and the Poisson error in comoving 
scale using the 1/U max method, where V maK is the volume 
corresponding to the total redshift range over which the 
galaxy would be detected at the K ~ 25 flux limit. 

4-5. Hosting Dark Halos and Halo Occupation Number 

Here we examine a characteristic mass of dark matter 
halos that hosting galaxies and the halo occupation num- 
ber from the galaxy clustering and number density com- 
pared with a theoretical model in the CDM cosmology. In 
Fig. 10 we plot the bias parameter b to compare the clus- 
tering strength with a dark matter halo model against the 
number density of the galaxies. In the previous section, we 
obtained the correlation length ro, fitting the power law 
ACF to the observations. The r value would be a rea- 
sonable quantity to compare the observations with dark 
matter models. However, it strongly depends on the as- 
sumed 7 and, moreover, could be affected by the small 
scale excess in 8 < 8" bins as shown in the previous sec- 
tion. For this reason, here we adopt the bias instead of the 
correlation length. In addition, for comparison, we depict 
the results of other studies based on if-selected catalogs 
(Daddi et al. 2003; Quadri ct al. 2007) and those of opti- 
cal samples (LBGs, Lee et al. 2006; DEEP2, Conroy et al. 



2007). 

The model predictions for the number density and aver- 
age bias of dark halos with mass larger than the minimum 
Mpg are shown by solid lines. The discrepancy in number 
density between the observed points and the model pre- 
diction indicates that the number densities of galaxies and 
dark halos are not exactly in one-to-one correspondence. 
We define the occupation number, N oc , as the ratio of the 
number density to that of the dark halos. 

Assuming that the bias of galaxies (b) reflects the bias 
of dark halo (&dh) hosting the galaxies (i.e., b ~ &dh), 
we estimate the hosting halo mass Mdh with the help of 
the analytic models given by Sheth, Tormen (1999) and 
Sheth, Mo, Tormen (2001) (e.g., see Ouchi et al. 2004 for 
the detail method) , where b is obtained by averaging the 
observed biases in the bins with 8" < 8 < 100" to avoid the 
small scale excess (see Table 1). For reference, we mark 
the Mdh values, which correspond to the bias, at the right 
side of Fig. 10 

We summarize in Table 2 the results of the minimum 
halo mass M^", the occupation number N oc , the halo 
mass -Mdh with the same bias as the galaxies (b ~ &dh), 
and the ratios of Mdh to the average stellar mass (M»). 

4-6. Redshift Evolution of Galaxy Biases 

In order to explore an ancestor-to-descendant connec- 
tion of different galaxy populations at various rcdshifts 
from the viewpoints of both stellar mass and dark matter 
assembly, we compare our results of the bias and those of 
literature with model predictions as a function of redshift 
in Fig. 11. If the bias values are not available in the liter- 
ature, we calculate the bias b g by the following equation 
from the published ro and 7, 

/ [8h-m P c/r }~ 
9 V &M(r = 8^Mpc) ! lUJ 

where £dm is the two-point correlation function of under- 
lying dark matter. We calculate ro of the dark matter us- 
ing nonlinear CDM correlation function (Peacock, Dodds 
1994). For comparison, we show the bias values of clus- 
ters of galaxies (Bahcall ct al. 2003) and groups of galax- 
ies (Padilla et al. 2004; Girardi et al. 2000) in the local 
Universe. The bias value of the local galaxies with typical 
luminosity M* ~ —20.5 (Zehavi et al. 2005), which corre- 
sponds to M* ~ 3.8 x 10 10 M Q using M*/Z ~3 (Kauffmann 
et al. 2003) , is also shown with the bias range for galaxies 
with -23 < M r < -17. 

If the motion of galaxies is purely caused by gravity and 
merging does not take place, the bias value of galaxies will 
decrease as the Universe evolves with time, according to 

b z = l + (b -l)/D(z), (12) 

where D(z) is the growth factor and 60 is the bias at 
z = (Fry 1996). We calculate the growth factor follow- 
ing Carroll, Press, Turner (1992). This galaxy-conserving 
model provides an upper limit of the bias evolution as 
shown in Fig. 11. 

On the other hand, if galaxies continue merging with 
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Table 1. Clustering measurements of galaxies with K < 25.0 
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Fig. 9. Same as Fig. 5, but for galaxies at 1 < z < 2 and 2 < z < 4, grouped with their stellar masses. 



the same merger rate of their hosting halos, the model 
gives a lower limit to the bias evolution. In order to com- 
pute the theoretical predictions for the bias evolution in a 
hierarchical structure formation model, we adopt the ex- 
tended Prcss-Shcchter formalism in the framework of the 
CDM cosmology with the power spectrum normalized to 
reproduce the local cluster abundance (as = 0.9) (e.g., see 
Hamana et al. 2006 for more details). Since the prob- 
ability distribution functions (PDFs) of the descendant 
mass are skewed towards larger mass, we adopt the mode, 
rather than the average, of the distribution as the typical 
mass and bias of the descendant halo (Af^ H and &p H in 
Table 2). It should be noted that the error of 6^ in- 
cludes a variety of descendant masses due to the different 
mass assembly history from the original halo mass Mdh 
at the observed epoch. The bias evolution paths of the 
typical halo hosting the galaxies of each subsample from 
the observed epoch to the present day are shown in Fig. 
11. Table 2 gives the bias error A&dh, which includes the 
errors of &dh and Mdh- The error bars (A&dh) shown in 
Fig. 11 arc considered to be the possible bias range of the 
descendants. 

5. Discussion 

We have measured the clustering properties of galax- 
ies in a 24.4 arcmin 2 area of GOODS-N region in rcdshift 
bins of 1 < z < 2 and 2 < z < 4, using our JT-selected 
catalog with the limiting magnitude K ~ 25 at the 90% 
completeness. Our results conclude that DRGs show the 



conspicuously large clustering amplitude at 2 < z < 4 and 
more luminous galaxies in K band are more strongly clus- 
tered. Small scale excesses (9 < 10", < 0.3 Mpc at z ~ 3) of 
the clustering amplitudes are discernible for the subsam- 
ples with 0.5 < J — K < 1.3, in which most of the LBGs 
populate. Since the flux in K band is well correlated with 
the stellar mass of galaxies at z < 4, the dependence of 
the clustering amplitude on K flux suggests that massive 
galaxies are more strongly clustered than low mass galax- 
ies. In galaxy evolution in the context of CDM models, 
the distribution of galaxies is strongly correlated with that 
of dark matter halos. The mass of the hosting dark halo 
could be one of the most fundamental quantities in galaxy 
evolution. Therefore, the evolution of the stellar mass of 
galaxies and dark matter should be investigated in a uni- 
fied way. In this regard, we divide our sample into two 
subsamples of massive (M* > 10 10 Mq) and less- massive 
(M* = 10 9 ~ 10 M Q ) galaxies, compare the bias with those 
of underlying dark matter and hosting dark halos. 

5.1. Descendants of the Galaxies 

We will discuss possible present-day descendants of the 
galaxies observed at 1 < z < 4 from the viewpoint of mass 
and the evolution of hosting dark halos. Comparing the 
clustering bias with the model predictions, we have es- 
timated the characteristic mass of hosting halos at the 
observed epoch (Table 2). Then, adopting the extended 
Press-Shechter model, we have computed the history of 
mass assembly of hosting halos to the present-day in terms 
of bias evolution (Fig. 11) to find out which populations in 
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Fig. 10. The relationship between galaxy number density and the bias of the galaxies at 1 < z < 2 (left) and 2 < z < 4 (right), 
divided into the groups by the stellar mass and color. The large open circles and squares shows the bias averaged in 8" < 9 < 100" 
bins for the massive (> 1O 1O M0) and low mass (10 9— 10 Mq) galaxies, respectively. The small filled and open circles (squares) depicts 
the color-selected subsample for massive (low mass) galaxies. For comparison, the results of FIRES (Daddi ct al. 2003), Quadri 
ct al. (2007), LBGs (Lee et al. 2006) and DEEP2 (Conroy et al. 2007) are depicted. The solid lines indicate the number density and 
average bias of the dark halo given by Sheth, Tormen (1999) integrated with the minimum mass (M™-^ n ) shown at the upper side of 
the figures with tick marks in units of Mq. At the right side for each figure, halo masses (Mqh) corresponding to b are marked. 



the local Universe are possible descendants. In particular, 
we will trace the evolutionary track of the predicted halo 
masses of DRGs, massive galaxies, and low mass galaxies 
at different epochs. Figure 11 suggests that the descen- 
dants of the galaxies possibly coalesce in halos of various 
masses, from those of normal galaxies to clusters of galax- 
ies. 

The evolutional track would reside in the range of 
galaxy-conserving model (upper limit) and merging model 
(lower limit) in Fig. 11. Taking account of the large er- 
ror budget of evolutional path of halos, we stress that the 
present-day descendants of DRGs are likely to be massive 
ellipticals located in groups of galaxies or clusters of galax- 
ies, This finding is in good agreement with the results in 
the previous clustering studies for DRGs (e.g., Grazian 
et al. 2006; Quadri ct al. 2007). In fact, the biases at 
various epochs arc roughly on the consistent evolutionary 
path (red-solid line in Fig. 11) from z ~ 3 to z ~ 1 within 
the error, except less massive DRGs sample of Grazian 
et al. (2006b) at z ~ 1.5. 

On the other hand, the bias values and the evolution 
for galaxies with M* > 1O 1O M and M* = 10 9 " 10 M o at 
1 < z < 2 and 2 < z < 4 trace the evolutionary paths of 
lower mass halos than those hosting DRGs. Although 
the error budgets on the bias estimate are very large, it is 
suggested that the galaxies with M* > 10 10 M Q evolve into 



luminous galaxies (< 3L*). LBGs are frequently argued as 
the strong candidates for the progenitor of the present-day 
early-type galaxies (e.g., Ouchi et al. 2004; Hamana et al. 
2006; Lee et al. 2006). The massive galaxies in the present 
sample would be a more general population for progenitors 
of luminous early-type galaxies. On the other hand, those 
with M, = 10 9 ~ 10 M Q are suggested to be the progenitors 
of less luminous local galaxies with M r ~ —20.0 (~ 0.6L*), 
which corresponds to M* - 2 x 10 10 M o with M*/Z ~ 2 
(Kauffmami 2003). In other words, the progenitors of 
present-day field galaxies were already populous at z ^ 
4 and bright enough to be observable with the limiting 
magnitude of K ~ 25. 

It should be noted, however, that only galaxies with 
M* = 10 9_10 M Q at 1 < z < 4 would not become present- 
day L* galaxies. The PDF of the present-day descen- 
dants has a very broad spread of mass range (Mf} H in 
Table 2) due to the wide variety of the mass assembly his- 
tory since observed epoch. Some different high-z galaxy 
populations may have evolved into a similar low-z popu- 
lation if the sample galaxies followed different evolutional 
paths. Since the PDFs are skewed towards larger mass 
(Hamana ct al. 2006), a certain fraction of lower mass 
galaxies (M* < 10 9 M Q ), which arc expected to be more 
numerous population at high-z, is likely to evolved into 
normal galaxies at the present-day. To learn what frac- 
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Fig. 11. The bias evolution with redshift for the sample galaxies in comparison with the previous observations at high- z, the results 
in the local Universe, and the CDM model predictions. The data in the relevant literature, which are represented in Figs. 6 and 8, 
are shown with the same symbols. In addition, optically selected samples for LBGs are referred to Lee et al. (2006) at z ~ 3.2 and 
Ouchi et al. (2004) at z ~ 4. The spectroscopic samples of DEEP2 are depicted at z ~ 1 (Coil et al. 2006). For the optical samples, 
the larger biases are for more luminous galaxies. The bias ranges for the samples at the local Universe (z ~ 0) are cited from Bahcall 
et al. (2003) for clusters of galaxies, Conroy et al. (2006) for groups of galaxies, Zehavi et al. (2005) for SDSS sample with the 
value for typical luminosity (M* ~ —20.5) in r band. The dotted lines show the bias evolution predicted by the galaxy-conserving 
model (Eq. 11). The merging evolution paths of the hosting halos in the framework of the CDM cosmology based on the extended 
Press— Shechter model are depicted by solid and dash lines for our samples. The error budgets of the bias evolution, which include 
the variety history of mass assembly from the observed epoch and the observational error of the bias, are depicted at the left lower 
corner with solid lines. 
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tion of lower mass galaxies contributes to the present-day 
normal galaxies, deeper observations in NIR will be nec- 
essary. 

5.2. History of Stellar Mass Assembly 

The dependence of dark-halo-mass to stellar-mass 
M^m/M* on galaxy properties and redshift is a clue to 
understanding the formation and evolution of galaxies in 
the context of merging and star formation in galaxies. 
The history of stellar mass assembly in galaxies have been 
studied by means of the stellar mass function (SMF) at 
various epoch in a statistical manner (e.g., Cole et al. 2001; 
Bell ct al. 2003; Rudnick ct al. 2003; Bundy ct al. 2006; 
Fontana et al. 2006). However, the SMF studies have been 
limited to comparatively massive galaxies at high-z. We 
present a schematic view of the stellar mass assembly his- 
tory since z ~ 4 to the present day from the viewpoint of 
the halo mass evolution of the galaxies. 

Conroy et al. (2006) have measured the virial-to-stcllar 
mass ratios for isolated ~ L* galaxies at z ~ 1 and z ~ us- 
ing the motions of satellite galaxies around isolated galax- 
ies in DEEP2 and SDSS catalogs. They gave M DM /M* ~ 
58-74 for z ~ galaxies (M* < lO n M ) (h = 0.7) with 
no strong dependence of the stellar mass. Provided the 
descendants of the galaxies observed in this study have 
Mrm/M* ~ 70 on average at z = 0, we obtain the av- 
erage stellar mass of the descendant using the (A/q H ) 
values in Table 2 to be (M*(jz = 0)) - 1.0 x 10 n M Q and 
- 2.3 x 10 n M Q for M* > 1O 1O M galaxies at 1 < z < 2 
and 2 < z < 4, respectively; the stellar mass increased 3-6 
times since observed epochs. Then the stellar accumula- 
tion rate, using the increased stellar mass, can be calcu- 
lated with the time interval between the observed epoch 
and the present day. The rates are 6.7M yr" 1 ((z) = 1.4 
to 0) and 17M yr" 1 ((z) = 2.8 to 0) for M* > 10 10 galax- 
ies. The SFR for massive galaxies at the local Universe 
is much smaller (Brinchmann ct al. 2004), so that general 
star formation would not account for the increased stellar 
mass. Some merger events would have occurred to the 
massive galaxies. 

In the same way, we evaluate the mass of descendants of 
10 9-io Mq g a i ax i es to be M„ ~ 1.6 x 10 10 M o , which is - 5 
times larger than that at the observed epoch. The stellar 
accumulation rates for the low mass galaxies are 1.4M 
yr" 1 ((z) = 1.4 to 0) and 1.1M yr" 1 ((z) = 2.7 to 0), 
respectively. The SFRs are even smaller than ~ 2M yr _1 
of the galaxies with M* ~ 10 10 A/ Q in the local Universe 
(e.g., Kewley et al. 2002; Brinchmann et al. 2004). If the 
SFRs continued constantly since z ~ 1 to the present day, 
the low mass galaxies are likely to have accumulated the 
stellar mass of the present-day normal galaxies without 
major mergers. 

It is interesting to note that the accumulation rate for 
massive galaxies is higher at 2 < z < 4 than at 1 < z < 2, 
while that of low mass galaxies tends to increase at lower 
redshift. In other words, the abrupt increase of stellar 
mass occurred during z = 2-4 by active star formation 
mainly in massive galaxies. This finding supports down- 
sizing of galaxy formation; more massive galaxies have 



higher SFR than low mass galaxies at high-z. 

Next, we examine major merger rate from the view- 
point of number density evolution of galaxies. The num- 
ber densities of M* > 1O 9 M galaxies are 6.6 x 10 _3 Mpc -3 
and 4.0 x 10~ 3 Mpc~ 3 at 1 < z < 2 and 2 < z < 4, re- 
spectively (Table 1). It increased 1.7 times during ~2 
Gy, which is also suggestive of active star formation in 
the era. The number density of the present-day galax- 
ies can be obtained using galaxy mass functions in the 
local Universe (e.g., Cole et al. 2001; Bell et al. 2003). 
The galaxies with M* > 1O 9 M at 1 < z < 4 evolved into 
M* > 5 x 1O 9 M galaxies in the present day as discussed 
above. Using the mass function of Bell et al. (2003) and 
Cole et al. (2001), we obtain the number density of the 
galaxies with M* > 5 x 1O 9 M at the local Universe to 
be ~8.8 x 10 _3 Mpc~ 3 . The similar or even larger density 
of the present-day galaxies indicates that major mergers 
are not significant events for the stellar mass accumula- 
tion. In the same way, we evaluate the merging rate of 
the massive galaxies (M» > lO lo Af ), which evolve into 
M» > 5 x 1O 1O M at the present day. The number den- 
sity of M* > 5 x 1O 1O M at the local Universe is (0.9- 
1.9) x 10~ 3 Mpc~ 3 , which depends on the mass functions 
of Cole et al. (2001) and Bell et al. (2003). It is compara- 
ble with the density 1.6 x lO^Mpc" 3 at 1 < z < 2 (Table 
1). Therefore, they arc likely to have experienced ~ 1 
major mergers since z ~ 1. 

However, the merger rate would be a lower limit if the 
lower mass galaxies (M* < 1O 9 M ), which were expected 
to be more numerous at high-z, evolved into the massive 
galaxies and give a significant contribution to the num- 
ber density in the local Universe. It is noted that star 
formation and merging arc degenerated in terms of stel- 
lar accumulation process. Further discussion would be 
beyond the scope of this study. 

5.3. Environmental Effect on Star Formation Activity 

We see the color dependence of the bias values in Fig. 
10. The U — V ~ 0.7 color corresponds to AO stars, 
so that the bluer color suggests active star formation. 
The red massive galaxies (M* > 1O 1O M ), which include 
DRGs, have a higher bias than that of the bluer sample 
at 2 < z < 4. (This is not confirmed at 1 < z < 2, because 
most massive galaxies are part of U— V ^ 0.7 populations.) 
However, the tendency is reversed for low mass galaxies 
at both epochs of 1 < z < 2 and 2 < z < 4; the blue low 
mass galaxies tend to have larger bias than the red galax- 
ies. The halo mass hosting low-mass bluc-galaxics and the 
occupation number (M DH , A oc ) arc (~3.5tf;g x 1O 12 M , 
^4.8^23) at 1 < z < 2, while those of the red low mass 
galaxies are (~l.l±f f x 1O U M , ~0-2±g;|). Also, at 
2 < z < 4, the bluer samples have larger host dark halo 
masses and occupation numbers, though the large error 
does not allow us to evaluate reliable values. The blue 
low mass populations are likely to coalesce in more mas- 
sive halos than red low-mass populations with higher halo 
occupation number. In other words, low mass galaxies 
in massive halos are more active in star formation than 
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those in lower mass halos; the fact is suggestive of the en- 
vironment effect by dark halo mass. The similar tendency 
has been observed for LBGs. The occupation number 
(A oc ~ 1) of LBGs in massive dark host halos (~10 13 M Q ) 
is much larger than that (A oc ~ 0.1) in lower mass dark 
halos (<lO 12 A/ ) (Ouchi ct al. 2004; Lee et al. 2006). 
The star formation rate of LBGs in massive halos with 
> 10 12 Mq are several to ten times larger than in less mas- 
sive halos (Cooray, Ouchi 2006). These facts also support 
the environment effect of dark halo mass on star formation 
activity. 
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